Constraining a spatially dependent rotation of the Cosmic Microwave Background 
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Following Kamionkowski (2008) , a quadratic estimator of the rotation of the plane of polarization 
of the CMB is constructed. This statistic can estimate a spatially varying rotation angle a(n). 
We use this estimator to quantify the prospects of detecting such a rotation field with forthcoming 
experiments. For PLANCK and CMBPol we find that the estimator containing the product of the E 
and B components of the polarization field is the most sensitive. The variance of this EB estimator, 
N(L) is roughly independent of the multipole L, and is only weakly dependent on the instrumental 
beam. For FWHM of the beam size Qf w hm ~ 5' — 50', and instrument noise A p ~ 5 — 50/xA'-arcmin, 
the scaling of variance N(L) can be fitted by a power law N(L) = 3.3 x 10~ 7 ApB}'^ hm deg 2 . For 
small instrumental noise A p < 5fj,K-a.rcmin, the lensing B-modes become important, saturating the 
variance to ~ 10 _6 deg 2 even for an ideal experiment. Upcoming experiments like PLANCK will 
be able to detect a power spectrum of the rotation angle, C aa (L), as small as 0.01 deg 2 , while 
futuristic experiment like CMBPol will be able to detect rotation angle power spectrum as small 
as 2.5 x 10~ 5 deg 2 . We discuss the implications of such constraints, both for the various physical 
effects that can rotate the polarization as photons travel from the last scattering surface as well as 
for constraints on instrumental systematics that can also lead to a spurious rotation signal. Rotation 
of the CMB polarization generates B-modes which will act as contamination for the primordial B- 
modes detection. We discuss an application of our estimator to de-rotate the CMB to increase the 
sensitivity for the primordial B-modes. 

PACS numbers: 



I. INTRODUCTION 



The polarization of the Cosmic Microwave Background 
(CMB) field can be studied in terms of the parity even 
E and parity odd B-modes [H, H, H, HI . In standard cos- 
mology, the physics governing the radiating field is par- 
ity invariant. Hence, the parity odd correlations (TB), 
(EB) vanish identically irrespective of the exact values 
of the cosmological parameters. However, the plane of 
linear polarization of CMB fields can be rotated due to 
interactions which introduce a different dispersion rela- 
tion for the left and right circularly polarized modes, 
during propagation from the surface of last scattering 
surface to the earth. Such rotations generate non-zero 
cross-correlations (TB), (EB) in the CMB field. Thus, 
measurement of these correlations allows us to estimate 
the rotation of the plane of the CMB polarization 
Such interactions can come from three main sources: (a) 
interaction with dust foregrounds, (b) Faraday rotation 
due to interaction with background magnetic fields, and 
(c) interactions with pseudoscalar fields [|| . The interac- 
tion with foregrounds leads to a a frequency dependent 
effect. The same is true of Faraday rotation, where the 
frequency dependence is (~ v~ 2 ) 0, H, H, E3] , while the 
interaction with pseudo-scalar fields is frequency inde- 
pendent. The distinct frequency dependencies allow one 
to separate these effects. 

We know that parity is violated by weak interactions, 
and is possibly violated in the early universe, to give rise 
to baryon asymmetry. Hence, investigating the existence 



of parity violating interactions involving cosmologically 
evolving scalar fields is well motivated. As an example we 
consider an interaction of the form Aj<' tl „F' u ' [flEl- It 
has been shown that such a term can rotate polarization 
vector of linearly polarized light by an angle of rotation 
a = jj J dr<j) during propagation for a conformal time 
t. The fluctuations in the scalar field 4> then will be im- 
printed in the rotation angle a of the polarization. It is 
interesting to ask what is the level of these fluctuations 
that can be detected with the upcoming CMB polariza- 
tion experiments. 

However, the observational situation is somewhat com- 
plicated by the fact that the measured CMB fields could 
be rotated with respect to the signal due to instrumen- 
tal systematics: a mis-calibration of the orientation of 
the instrument which results in a constant rotation and 
differential offsets of the orientation of the individual de- 
tectors of the instrument resulting in rotation dependent 
on angular position n. It should be noted that this sys- 
tematic effect is also a concern for the detection of po- 
larization B modes, a major goal for forthcoming polar- 
ization experiments, since it can result in a spurious B- 
mode detection. Rotation of the CMB a(n), either due 
to interactions with a pseudoscalar field, or due to the 
instrumental rotation miscalibration can be a function of 
angular position (n). 

An estimator for the spatially varying ration angle 
a(n) has been reported in Kamionkowski [l2l |. Here, 
we point out that the widely used formalism for gravita- 
tional lensing [k| may be suitably modified to describe 
the rotation of CMB polarization. We use this to con- 
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struct approximate, but simpler form of the quadratic 
estimator of a(n) using the flat-sky limit, and study its 
variance. 

These estimators may be used to study the physics 
behind the rotation of the polarization of light. As we 
shall discuss, we can also put an upper-bound on the fre- 
quency independent rotation from non-standard interac- 
tions. We also discuss the use our the estimator d(n) 
to control instrumental rotation systematics for the de- 
tection of primordial B-modes. The presence of rotation 
systematics in the instrument generates B-modes. Hence, 
the control of rotation systematics of instruments is im- 
portant for the measurement of B-modes of the CMB 
polarization; a major goal of subsequent polarization ex- 
periments. In this paper, we will discuss the prospects 
of detecting non-standard physics through measuring the 
angle of rotation, and the level to which we can control 
rotation systematics by this estimator. 

A search for a constant rotation of the polarized light 
by an angle a from radio gala xies and the CMB is already 
underway 0, EI EE EMl EE HE EH ■ So far, there is 
no evidence of non-zero angle of rotation, and the angle 
a is constrained to be less than a few degrees ED, EH • At 
present, there are no studies of constraints on spatially 
varying rotation angle a(n). 



II. FORMALISM 

In this section, we will construct an estimator for the 
spatially varying rotation field. We shall also describe a 
physical scenario which give rise to a frequency indepen- 
dent but spatially varying rotation. 



Estimator for Spatially Dependent Rotation 
Field 



We will describe the observable effect of rotation on 
the CMB polarization fields. Let the un-rotated (usual) 
CMB temperature field and the Stokes parameters at an- 
gular position n be T(n), and Q(n), U(h) respectively 
The relevant ensemble averages of the un-rotated CMB 
field can be encapsulated in 

= 0, (x*(l)x'(l')) = (2nf5(l - l')Cr' , (1) 

where x, x' run over the T, E, or B fields, and Cf x is 
the un-rotated CMB power spectrum. The temperature 
fields are invariant under a rotation of the polarization 
by an angle a (A) at the angular position n, while the 
Stokes parameters transform like a spin two field. Thus, 
due to rotation, the observed fields are 

(Q(n) ± iU(n)) = (Q(n) ± U(h)) exp(±2m(n)). (2) 

The E and B fields of the CMB can be constructed from 
observed Stokes parameters. In a Fourier basis, in the 



TABLE I: Minimum variance filters 

xx' f xx '{h,h) 

TE 20?^ sin2<£ lll2 

TB 2<5£ B cos2y>ni 2 

EE 2[C^-C i f B ]sin2 Wll2 

EB 2[C ££ -C i flfl ]cos2^ lll2 

BB [C i f B +C« B ]sin2^ ll2 



flat sky approximation, 



[E±iB){\) = dh[Q(h)±iU(h)] 



,^2iipi — iln 



(3) 



where ip\ = cos~ 1 (n • 1). Since, the angle of rotation is 
already constrained to be small, we will work out the 
effects to first order in a(n). 

Since, we can only compute correlations of the CMB 
polarization modes theoretically, we want to isolate the 
change in correlations due to this rotation. Even in the 
absence of the physics causing rotation, we expect the 
CMB fields to be gravitationally lenscd by matter in- 
homogeneities. Hence, the change due to rotation is 
the difference between lensed rotated fields, and lensed 
un-rotated fields T, E,B. We make the similarity of 
our problem with gravitational lensing of the CMB [l3| 
manifest by writing the change in the CMB field modes 
5x(l) = x(l) — x{V) due to rotation 



0. 



ST(l) 
SB(l) 

SE(l) = - 



fill r 



(27T) 2 

d 2 V 



(4) 

E(Y) cos 2ip n - B(l') sin 2</?iJ W(L), 
B(l') cos 2ip V i + E{V) sin 2tp vl ] W(L), 



where L = 1 — 1', and ipw = (pi — tpy and 1T(L) = 2a(L). 
Thus, due to rotation, a mode of wavevector L mixes 
the polarization modes of wavevectors 1 and 1' = 1 — L. 
Taking the ensemble average of the CMB fields for the 
fixed a field, one gets 



(^(l)x'(l')) 



CMB 



(x*(l)£'(l')>+/xx'(M'HL). (5) 



The TE, and TB correlations are produced indirectly via 
non-zero primordial TE correlation. 

Our goal is to use Eq. ([5]) to construct a suitable 
estimator of the Fourier components a(L) of the rota- 
tion field in terms of the observed fields T(l), E(l), B(l) 
and a theoretical computation of the power spectra in- 
volving un-rotated fields. Following Hu and Okamoto 
EH, we can define an unbiased quadratic estimator 
a xx '(L) for a(L) for each combination of the CMB modes 
x(li),x'(l2) by weighting quadratic combinations of dif- 
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ferent polarization modes by F xx i (li, I2) appropriately: 

a ra ,(L) = A XX ,(L) J ^-lx{\ x )x'{h) 

-(x(h)x'(h)}} F xx ,(h,l 2 ), (6) 
where L = I2 — li, and the normalization 



A xx ' (L) 



d 2 h 



-I -1 



(7) 



is chosen to make the estimator unbiased, i.e. (d(L)) = 
a(L). The fields x(l) can be obtained from the map of 
an experiment, while the CMB power spectrum of un- 
rotated but lensed fields can be computed from publicly 
available Boltzmann codes like CMBfast and CAMB. 

The weights F xx i can be optimized by minimizing the 
variance (a xx >(L)a* x , (L')). For xx' = TB and EB 



Fxx'(h, b) 



fxx'{h,h) 

qxx^jx'x' ' 



(8) 



where Cf* and Cf 2 x are the observed power spectra in- 
cluding the effects of both the signal and the noise, 



Ci 



(9) 



where C% x is the noise power spectrum. We assume 
the detector noise to be known apriori, be isotropic 
and Gaussian distributed. We include effects of beam- 
smearing by a symmetric Gaussian beam. Then, the 
noise power spectrum is 



= A 



J8\n2 



(10) 



where A x is the instrument noise for temperature (x=T) 
or polarization (x=P); and Qfwhm is the full-width half- 
maximum (FWHM) resolution of the Gaussian beam. 
We will assume a fully polarized detector, for which 
V2A T = A P . 

The variance of the estimator is 

(a xx ,(L)a* xx ,(L>)) = (2nf5(L - L'){Cg« + N XX ,(L)} , 

(11) 

where for the minimum variance estimator, the Gaussian 
noise N XX /(L) = A xx r(L), and gives the dominant contri- 
bution to the variance. 

The Gaussian noise N(L) is dependent only on the 
instrumental noise power spectrum Cf x ™ and the power 
spectrum of the un-rotatcd polarization field Cf x . Hence 
the estimator noise depends on the cosmological param- 
eters through the power spectrum of the polarization 
fields. We choose a standard fiducial model with a flat 
AC DM cosmology, with no rotation (i.e. a = 0), with 
parameters described by the best fit to WMAP5 [Til ], 
given by fl b = 0.045, fl c = 0.23, H = 70.5, n s = 
0.96, n t = 0.0, and r = 0.08. 



Since, in reality the polarization field will be gravita- 
tionally lensed by the inhomogeneities in the matter dis- 
tribution in the fiducial model, it is appropriate to use the 
power spectrum of lensed anisotropies as the un-rotated 
field in calculating N(L). An angular remapping of pho- 
ton positions due to gravitational lensing may result in 
an apparent frequency independent rotation of the plane 
of polarization, potentially biasing the estimator. Here, 
we show that the bias is negligible. Taking lensing into 
account the average of the estimator is 

{a xx '(L))cMB = 
a(L)+A xx ,(L) J ^Lf^F xx ,(h,h)<t> len (L) 7 (12) 

where cf> len is the line of sight projection of the gravita- 
tional potential ^(x). The first term on the right hand 
side is the desired rotation field, and the second term 
represents the bias from lensing. The form of lensing fil- 
ters f l x e £ s can be found in Table I of [TJ[. The lensing 
filter f x e ™ are nearly orthogonal to the rotation window 
F xx /. Hence the integrand of the lensing bias oscillates 
around zero, and even for cf) len ~ 1, the bias is negligi- 
bly small in comparison to the square root of Gaussian 
noise \J N(L). Since Gaussian noise sets the minimum 
detectable rotation a(L), we can neglect the bias for all 
practical purposes. 



B. Rotation from non-standard interactions 

As discussed in introduction, apart from instrumen- 
tal systematics, which we will discuss in section UlI Bl a 
way of generating frequency independent rotation is by 
pseudoscalar fields coupled to photons. 

There are no pseudoscalars in the standard model of 
particle physics that couple to radiation. However, they 
are common in particle physics beyond the standard 
model [H, H3, 0, HH- In cosmology they have been 
invoked for dark matter or dark energy models, and also 
as solution to the fine-tuning problem of dark energy 

A pseudoscalar field <fr can couple to the electromag- 
netic fields by a Chern-Simons interaction term [ill . 132 . 




C 



1 



1 



2M 



(13) 



where F^ v is the electromagnetic field strength tensor, 
F^ v is its dual; and the pseudoscalar coupling with elec- 
tromagnetic field is supressed by mass scale M . 

The interaction term ^jF^F' 11 ' is invariant under 
U(l) gauge transformations, Lorentz symmetry and par- 
ity, and is suppressed by a mass scale M. The fact that 
no such effect has been detected in the laboratories puts 
a lower bound on M. It has been shown that such a term 
can rotate polarization vector of linearly polarized light 
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by an angle of rotation a = -h J drcj) during propaga- 
tion for a conformal time r @ , which is largest when the 
source of polarization is farthest, which happens for the 
CMB @. 

The angle of rotation along a line of sight depends on 
the change of </> along that line of sight. We can write 
the field </> in terms of a homogeneous piece 4'q{ t ) an( i 
a position dependent perturbation S(f>(h, r). Then, the 
rotation angle of the CMB polarization 

a ^ = = M^ 0< - r °^ ~ ^o(Tdec) - S(j)(n, T dec )} , 

(14) 

where r and T dec are the conformal times today and at 
the surface of last scattering respectively, and the per- 
turbation S(f> today, at the detector can be taken to be 
zero without any loss of generality The shift symmetry 
of the Lagrangian implies that the field (f> is classically 
massless. In our toy example we will assume that this 
was the case during inflation, and the quantum fluctu- 
ations during inflation were frozen in the field and will 
result in spatial fluctuations today. In this case we can 
write the power spectrum of fluctuations as nearly nearly 
scale invariant, 

<0(k)^(k')) = (27r) 3 P(fc)<5(k-k') 

= (2 7 r) 3 c <?i /c™ 3 " 4 ( 5(k-k') (15) 

with n s = 0.96 [3, and = H 2 /2, where H is the 
Hubble parameter during the inflation. Computing the 
transfer function we find the power spectrum of the ro- 
tation angles to be 

Cr = J^J k 2 dkP^k)j 2 (kr)A 2 (k, T dec ) , (16) 

where r = To — T dec is the distance to the surface of 
last scattering, and A^(fc, r^ee) = 3ji(kT rec )/kTd ec is the 
transfer function which is unity for scales larger than 
horizon size and decay with an oscillating envelop for 
modes inside the horizon. 

We have discussed the scenario of physical interactions 
leading to rotation of the polarization field. There are 
two important physical parameters which determine the 
magnitude of the effect in this toy model: (a) the energy 
scale (or equivalently the Hubble rate) during inflation, 
which sets the amplitude of the fluctuation power spec- 
trum of mass-less fields, and (b) The mass scale M by 
which the Chern-Simons interaction term is suppressed. 
The mass scale M is certainly much larger than the cur- 
rent energy scales of particle physics experiments ~ 10 2 
GeV, and could be around or higher than the GUT scale 
~ 10 16 GeV. In fact, the current best constraints on 
mass scale M come from the upper limit on B-modes, 
M > 2.4 x 10 14 iJi4 GeV with H 1A = if/10 14 0. The 
exact value of the Hubble scale during inflation is also un- 
known, with current estimates suggesting it to be around 
~ 10 14 — 10 15 GeV. For reference the tensor to scalar ratio 



of perturbations produced in standard models of inflation 
is given by r = 0.14 H 2 A . This for current WMAP5 con- 
straints [Tj] of r < 0.2 translates to Hu < 1.2. Since 
the potential of detection is linked to an assumed energy 
scale of inflation, and M, it is useful to write rotation 
angle in terms of H/AI. The rms value of rotation angle 

for L 100 can written as J L c ~ IM. deg. Our 
approach is to find the magnitude of the rotation which 
would be detectable. 



III. RESULTS 
A. Estimator 

We discuss the prospects of using our unbiased esti- 
mator d(n) for detecting Fourier modes of the spatially 
varying rotation a(n). We study the estimator variance 
and its dependence on instrumental characteristics focus- 
ing on planned missions. 

We want to find out the magnitude of rotation that 
would be detectable. For this, the rotation signal must be 
larger than the noise in the estimator N xx > (L). In Fig. [IJ 
we show the variance for the EB (red curves) and TB 
(blue curves) estimators 1 as a function of multipole L 
for the three experimental setups, (1) PLANCK satellite 
with noise level Ap = 56^X-arcmin and FWHM of 7', 
(2) experiment with noise level Ap = 9.6/iif-arcmin and 
FWHM of 8', typical of upcoming ground and balloon- 
based CMB experiments (hereafter called Expl), and (3) 
a CMBPol-like instrument with noise level Ap = \[2[iK- 
arcmin and FWHM of 4', typical of future space-based 
CMB experiments. The black solid line shows the power 
spectrum of the cosmological rotation signal Cf a from 
the Chern-Simon coupling discussed in section III Bl We 
show this signal for c^/M 2 = 10~ 4 , which corresponds to 
71/ = 1.5 x 10 15 i?i4 GeV. 

For the three experiments in consideration, the EB 
estimator is found to be the most sensitive, and the vari- 
ance is roughly constant up to L ~ 1000. Although, we 
do not show L = in the plot, our estimator can also 
be used to estimate the dctectability of uniform rotation. 
For uniform rotation, L = li — I2 = in Eq. (J7J), hence 
there is no mode mixing between different wavevectors. 
The variance N(L) for monopole and dipole are com- 
parable; so rotation a(L) can be constrained to similar 
levels for these modes. 

Currently only the constant angle of rotation (i.e. 
L = of our estimator) has been constrained, a < 2° 
ITU : and at present there arc no constraints on the 
spatial variation of the rotation angle a(n). From our 
Fig.Q]upcoming experiment like PLANCK will be able to 



1 For the three experiments considered in Fig. [T] the other estima- 
tors have much smaller sensitivity to the rotation angle. 
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FIG. 1: Estimator variance N(L) for the EB (red curves) and TB (blue curves) estimator as a function of multipole L. We 
show the noise for three experimental setups, CMBPol (dot-dashed), Expl (dot) and PLANCK (dashed). The solid black 
curve shows the rotation angle power spectrum for the model in which the pseudoscalar coupling to the electromagnetic field is 
suppressed by a mass scale M [Tl|], and the perturbations in <j> are seeded during the inflationary phase. We choose the fiducial 
value of the amplitude for power spectrum c^/M 2 = 10 -4 . This correspond to energy scale M ~ 10 15 -Hi 4 GeV, where -Hi 4 is 
the Hubble parameter during inflation in units of 10 14 GeV. Note that with CMBPol like experiment, one is sensitive to energy 
scale as large as M ~ 10 17 -H14 GeV, i.e. three orders of magnitude of improvement over the current best constraints. 



detect the rotation angle power spectrum Cf a as small 
as 0.01 deg 2 , while futuristic experiment like CMBPol 
will be able to detect rotation angle power spectrum as 
small as 2.5 x 10~ 5 deg 2 . These numbers translate to 
minimum detectible H/M = 4 x 10~ 3 for PLANCK, and 
H/M = 2 x 10~ 4 for CMBPol. 

In Fig. [21 we study the dependence of the estimator 
variance with the detector noise A, and the beam size 0. 
As in Fig. [TJ here also we show the cosmological signal 
C aa (L) for reference, with the same fiducial parameter 
c^/M 2 = 10~ 4 . The upper panel shows N EBtEB (L) as 
a function of L for various choices of instrumental noise 
A p but for a fixed beam size of 4'. The lower panel 
shows estimator noise Neb,eb(L) as a function of L for 
A p = 2fiK-axcmm and for various choices of beam size, 
starting from 20' to 1'. For FWHM of the beam size 
Qfwhm ~ 5'— 50', and instrumental noise A p ~ 5— 50/iif- 
arcmin, the scaling of variance Neb,eb(L) can be fitted 
by a simple power law N E b,eb{L) = 3.3 x 10 _7 A 2 9 1 c -^ m 
deg 2 . However, we cannot minimize the estimator noise 
N(L) to arbitrarily low levels by reducing the detector 
noise; the estimator variance plateaus out at a level of 
A p < 5/i-ftT-arcmin (larger than detector noise levels in 
CMBPol) to w 10~ 6 deg 2 . Physically, this is due to lens- 
ing effects. 

As discussed in Sec. Ill A[ the estimator variance 
N xx ',xx'(L) shown in Fig. [T] includes contribution due 
to lensing of E modes to B-modes. For experiments 
like Expl and PLANCK with A p > 5/iK-arcmin, lens- 



ing effects are negligible compared to the detector noise. 
Therefore, there is no difference in the estimator variance 
if the lensed polarization fields used as the un-rotated 
fields in Fig. [TJ are replaced with un-lensed polariza- 
tion fields. On the other hand, for experiments with 
A p < 5/iK-arcmin (like CMBPol), lensing of the CMB 
power spectrum dominates the estimator variance, and 
eventually limits the sensitivity of the an idealized instru- 
ment to ~ 10~ 6 dcg 2 . Further, wc calculated the leading 
order lensing contribution to noise, N lens -(L) which is 
proportional to lensing power spectrum C len (L), and is 
related to the connected part of the trispectrum. We find 
that this noise N lens -(L) is smaller than the estimator 
noise N xx i >XX '{L) shown in Fig. [T] for all cases. 

Can Estimator noise be further reduced? Lensing B- 
modes can be measured and hence, in principle, can be 
separated (de-lensing) from the pure rotation a(h) con- 
sidered above. In the absence of lensing, the sensitivity of 
the idealized instrument would be limited by the cosmic 
variance of primordial B-modes. While de-lensing can 
improve the sensitivity of the idealized instrument, to the 
level of the noise inherent in the de-lensing process, it is 
likely to be challenging [33|. In Fig. [3] we show how much 
the Gaussian noise for the EB estimator N(L)eb,eb re- 
duces as a function of amount of de-lensing. The de- 
lensed B-modes 

j^de—lens j j^lens (17) 

are used in the estimator which reduces the variance, de- 
pending on the fraction of de- lensing /. The lower curve 
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FIG. 2: Dependence of the variance N(L) of the EB estima- 
tor on instrumental characteristics as a function of multipole 
L: Upper panel: The diagonal lines represent the variance 
N(L) for fixed FWHM of 4' but varying detector moise A p . 
Lower panel: The diagonal lines represent the variance N(L) 
for fixed detector noise of A p = 2fj,K-a,rcmin but varying the 
beam size. In both the panels the oscillatory curve repre- 
sents the power spectrum of rotational field for the fiducial 
amplitude of c,p/M 2 = 10 -4 . 



(solid red) shows the most optimistic scenario where all 
the CMB modes are being de-lensed (up to L = 3000). 
The upper curve shows, although still challenging, a more 
realistic case where only the modes L < 500 are being dc- 
lenscd. 



B. Detectability of cosmological rotation and 
instrumental systematic effects 

As indicated before, a cosmological rotation field cv(n) 
can be confused with the instrumental rotation system- 
atics. A calibration error in the angular position of the 
instrument is degenerate with a spatially constant an- 
gle of rotation (L=0 of a(L)), while errors in the rota- 



FIG. 3: Gaussian noise (at L=2) for the EB estimator as 
a function of fraction of De-lensing of the CMB. Since the 
lensing effects are only important for experiments with A p < 
5/iK-arcmin, we assume an idealized experiment (i.e. A = 
0, Q fmhm = 0) to see the effect of de-lensing. Lower (solid red) 
curve assumes that all the observed B-modes (l m ax = 3000) 
are being de-lensed by equal amount. Upper (dash black) 
curve shows the effect of de-lensing when only the modes with 
£ < 500 are being de-lensed, while no de-lensing for I > 500. 



tion calibration of individual detectors in the instrument, 
leading to a relative mis-alignment of axes of the individ- 
ual detectors by angles tUi are degenerate with spatially 
varying a(n). For an instrument with a large number 
of detectors, we can treat the angles of rotation of the 
i th detector u>i as a smooth field as a function of the de- 
tector position. The relative offsets in the polarimctcrs 
in the detector could result in a systematics signal u>(n) 
in the map if the weighting of each polarimctcr changes 
from pixel to pixels in the map. This depends on the 
scan strategy. For illustration purposes we can model 
the statistical properties of systematics signal [HI, HH as 
a statistically isotropic Gaussian field with a power spec- 
trum given by 

C uu = cxp(-/(f + l)^) (18) 

where characterizes the rms value of this field w, and 
(T(jj is a coherence length. 

To see the effect of rotation systematic field, we can 
change a(n) in Eq. Q to a(n) + w(n), and re-dcrivc 
our estimator. The systematics field biases our estimator 
d(L) by an amount w(L) and increases the variance of 
the xx' estimator by an amount C WW (L) + N^, XX ,{L) 
(see appendix). For our model of systematics field and 
assuming that a and u> fields are uncorrelated, the power 
spectrum of bias is given by Eq. (jig)) i.e. C w "(i). For 
the expression for the systematic noise eb(^) pl casc 
refer to the appendix. In order to use the estimator for 
detection of a rotation field with power spectrum C aa (L) , 
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FIG. 4: Left Panel: Dot dashed (blue) curves show the systematics rotation power spectrum C UUJ (L) for three combinations of 
rms amplitude A u (in arcmin) and FWHM of coherence length, -^8 ln(2)<r UJ (in arcmin). The two dashed red curves show the 
Estimator variance Neb,eb(L) for the PLANCK (upper) and CMBPol (lower) experiment. The black oscillatory curve shows 
the power spectrum of rotational field for the fiducial amplitude of c^/M 2 = 10~ 2 , 1CP 4 , and 10~ 6 . Right Panel: For CMBPol 
experiment, solid red curve shows the estimator variance Neb,eb{L) and the dot dashed (blue) curves show the systematics 
noise N%% tEB (L) for three combinations of rms amplitude A u (in arcmin) and FWHM of coherence length, y/8 ln(2)er w (in 
arcmin) as the left panel. 



C aa (L) » C U "{L). 

In order for the cosmological rotation field a(n) to be 
determined to the noise levels in Fig. [I] both the system- 
atic bias power spectrum C^(L) and systematic noise 
^xx> xx' W should be smaller than the estimator noise 
N xx ',xx'{L). In Fig. [4] we study the dependence of the 
systematic field power spectrum C UUJ (L) and the system- 
atics noise term eb(L)- 

In the left panel of Fig. [4] the dot dashed (blue) curves 
show the systematics rotation signal for various choices 
of rms a mplitude and FWHM of coherence length 
y/8 ln(2)er. The systematics signal does not depend on 
the instrumental noise A, and beam, Qf w hm- Solid 
(black) curves show the fiducial cosmological rotation 
power spectrum C aa (L) for C(j> /M 2 = 1CT 2 , 1CT 4 , &1CT 6 , 
and the two dashed (red) curves show the Gaussian noise 
for the EB estimator N E b,eb{L) for PLANCK (upper) 
and CMBPol (lower) experiments. 

In right panel of Fig. 2] the dot dashed (blue) curves 
show the systematics noise for EB estimator eb(^) 
for the CMBPol like experiment and for various choices 
of rms amplitude A u and coherence length a. The solid 
(red) line shows the Gaussian noise for the EB estimator 
Neb,eb{L) for CMBPol experiment. 

Note that like the Gaussian noise Neb,eb{L), the sys- 
tematic noise N^g eb(^) 1S a ^ so weakly dependent on L. 
However the systematic power spectrum (depending on 
the coherence length) may only be similarly flat only to 
about L ~ 100 (for coherence 60'). Hence, we may char- 
acterize their values by their values at a particular value 
of L < 100. Together, the two panels show that at a co- 
herence length of 10', an rms amplitude of w 0.01' gives 



a systematic variance which is about ten times smaller 
the estimator variance in the CMBPol experiment, while 
the bias power spectrum is only about half the estima- 
tor variance. Smaller coherence length a and rms ampli- 
tude of the systematic fields result in smaller effect of the 
rotation systematics. This implies, that in order to de- 
tect rotation, using the CMBPol experiment, one would 
have to control the systematic field to much better than 
{cr w , A,} values {10', 0.01'}. 



C. De-rotating CMB to improve the sensitivity for 
the Primordial B-modes Detection 

An important design goal of the futuristic CMB po- 
larization experiments is the detection of primordial B- 
modes. Rotation generates B-modes (via Eq. (4)) which 
can be confused with the primordial B-modes. Both the 
instrumental rotation systematics and any cosmological 
rotation will generate B-modes. Hence, in order to study 
primordial B-modes, it is necessary to know the level of 
these spurious B-modes. 

An important application of our estimator is to mea- 
sure rotation and then in turn de-rotate the CMB polar- 
ization field to remove the spurious B-modes and hence 
increase sensitivity to the primordial B-modes detection. 
For this application, it is not important to know what 
the source of this rotation is. 

A specific example is the case when cosmological ro- 
tation is known (or assumed) to be small and we are in- 
terested in controlling instrumental systematics to detect 
primordial B-modes. One can in this case layout specifics 
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FIG. 5: The requirement on control of the systematic fields w for CMBPol {Left Panel), and Expl [Right Panel) for the 
detection of primordial B-modes, and the level to which the CMB can be de-rotated.. The dashed black line shows the 
systematic parameters {A, a} which generate spurious B-modes of the same magnitude at (L = 90) as the primordial B-modes 
(at L=90, where the primordial B-modes peak); for values of r = 0.1 (upper back) and r = 0.01 (lower black). The solid red 
curve represents the systematic fields for which the variance contribution due to the systematic field Neb.eb(L — 2) becomes 
equal to the estimator variance Neb,eb(L = 2). The blue line represents the systematic fields for which the power spectrum 
of the ui field at L = 2 is equal to the estimator variance Neb,eb(L = 2). 



on what is the minimum B-modes amplitude that will 
be detectable without worrying about rotation system- 
atics. The amplitude of the primordial B-modes is fixed 
by amplitude of tensor perturbations which depends on 
the energy scale of inflation. Equivalcntly one can use 
the ratio of amplitude of tensor and scalar perturbations 
r to characterize the B-modes. 

In Fig. [5] we show the comparison of required control 
of the systematics ui for CMBPol (left panel), and Expl 
(right panel) for the detection of primordial .B-modes, 
and the level to which the CMB can be de-rotated. The 
dashed black lines show the A u , and a u of the systematic 
fields (see Eq. (|18p ) for which spurious £?-mode power 
spectrum C uui {L) at L = 90 is equal the primordial B- 
modes power spectrum (at L=90) for values of r = 0.1 
(upper curve), and r = 0.01 (lover curve). The solid 
red curve shows the A u , and <r u of the systematic field 
for which the variance contribution due to the systematic 
field Neb EB (L) at L = 2 becomes equal to the Gaussian 
noise Neb,eb(L) at L = 2. The dot-dashed blue curve 
shows the A^, and a w of the systematic field for which 
the power spectrum of the u> field C UU (L) at L = 2 is 
equal to the estimator variance N^g eb(^) a ^ L=2. 

In order to be able to detect £>-modes of a particular r 
value, the parameters {A, a} must be in the region below 
the corresponding black lines. Dot-dashed blue curves 
represent the level of rotation systematics to which our 
estimator can be used to correct for; i.e. if the systemat- 
ics, {A, a} are above the blue curve then our estimator 
can reduce it to the level given by blue curve. If the in- 
strumental systematics parameter {A, a} arc above the 
red curve, our estimator is no longer a minimum variance 



estimator, and one should include the systematics noise 
effect N uu} in the variance of our estimator. However if 
instrumental systematics parameter {A, a} are below the 
red curve, then for a given r if the systematics require- 
ment for B-mode is less stringent that the requirement 
from rotation, one can use our estimator to control the 
rotation systematics for B-modes detection. 

If the systematics parameters {A, a} are below the dot- 
dashed blue curve, then the rotation from the systematics 
is smaller than the sensitivity of the estimator. If our 
estimator detects rotation signal in this case then it can 
be attributed to cosmological signal or to incorrectness 
of the systematics model. 



IV. SUMMARY AND DISCUSSION 

There are interesting physical mechanisms that can ro- 
tate the plane of CMB polarization. We presented ex- 
plicit formulae for estimators of the spatially varying ro- 
tation angle a(n) that can be constructed from future 
datasets in the flat sky limit. By computing the variance 
of these estimators, we estimate how large a variation 
in the angle a must be to be detected by a particular 
experiment. Currently only the constant angle of rota- 
tion (i.e. L = of our estimator) has been constrained, 
a < 2° fl4l [l5j. At present there are no constraints on 
the spatial variation of the rotation angle a(n). From our 
Fig.rfjupcoming experiment like PLANCK will be able to 
detect the rotation angle power spectrum Cf a as small 
as 0.01 deg 2 , while futuristic experiment like CMBPol 
will be able to detect rotation angle power spectrum as 
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small as 2.5 x 10 5 deg 2 . These numbers translate to 
minimum detectible H/M = 4x 10~ 3 for PLANCK, and 
H/M = 2 x 10" 4 for CMBPol. 

Gravitational lcnsing does not bias the estimator, how- 
ever it increase the variance of the estimator. The in- 
crease in the variance is sub-dominant for experiments 
with A p > 5/iK-arcmin. For small instrumental noise 
A p < 5/^/Garcmin, the lcnsing B-modes become impor- 
tant, saturating the variance to ~ 10 _6 deg 2 even for an 
ideal experiment. 

The physical mechanisms that give rise to the rotation 
field all probe interesting cosmological physics; in princi- 
ple, they can be separated using their frequency depen- 
dence and used to study the magnetic field or polariza- 
tion dust maps. A cosmological source for a frequency 
independent pure rotation field can be interpreted as ex- 
otic feature signifying a clear departure from standard 
model physics. Such a departure could be a violation of 
the equivalence principle, or violation of Lorentz invari- 
ance |36l . |37| . In the context of cosmology, an important 
example is pseudoscalar fields which have been proposed 
as dynamical models of dark energy to solve the cos- 
mological constant problem, dark matter, and also as a 
solution to the fine tuning problem of dark energy. 

We have considered a model in which the perturbations 
in the scalar field were imprinted during inflation and the 
scalar field couples to photon via Chern-Simons coupling 
discussed in Sec. Ill Bt which are suppressed by mass scale 
M. With CMBPol like experiment one can constrain the 
mass scale M > 10 17 H U GeV, where H 14 is the Hubble 
parameter during inflation in the units of 10 14 GeV. 

If the scalar field is assumed to be responsible for the 
dark energy, the constant rotation would probe the cou- 
pling scales of such a dark energy field, and establish its 
dynamic nature. The spatially varying part of this rota- 
tion field would also probe the clustering of such a field. 
Typically, such fields would have large sound speeds, so 
that clustering is only possible at large scales. Hence, 
the possibility of detection of the spatially varying field 
is best at low multipolcs. If one assumes that the field 



has the kind of clustering discussed in (38|, then from 
an experiment like CMBPol, one can constrain the mass 
scale of suppression of the Chern-Simons coupling term, 
to M > 10 10 GeV. 

However, the rotation field induced by a cosmological 
source can be degenerate with the rotation systematics 
of the instrument, which arc limited by the rotation cal- 
ibration of the polarimcters. Thus, the detection of such 
cosmological signals is only possible if the cosmological 
signal is larger than the level of rotation systematics sig- 
nal that can be controlled. Rotation systematics also 
effect the variance of our estimator. We quantify the 
level of systematics control required for detection of the 
cosmological signal to be only limited by the estimator 
noise N(L). If one can show from other experiments that 
the sources of such cosmological signals can be limited to 
magnitudes |<5a(L)| smaller than these systematic levels, 
then we cannot detect the cosmological signals. However, 
then we can use this fact that the observed rotation field 
should be less than this magnitude to calibrate the in- 
strument to control the level of rotation systematics to 
precision levels of ~ |<5a(L)|. This could enable a better 
study of effects like primordial B-modes, lensing, or the 
frequency dependent signals from sources like magnetic 
fields or foreground dust. Thus precise studies of this 
rotation field could either probe exciting physical effects, 
and/or enable better control of calibration and instru- 
mental statistics. 
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APPENDIX A: CONTAMINATION 



Instrumental rotation systematics and lensing of the 
CMB can effect our estimator. In this appendix we will 
show how the rotation systematics and lensing of CMB 
appear in our estimator. Let us denote systematics field 
by lo(L) and lensing filed by 4> len (L). We can incorporate 
the effect of rotation systematics by changing a(n) to 
a(n) + w(n), and effect of lensing by changing n to n + 
d(h) in Eq. (1), where d = V</>' e ™. In the presence 
of rotation systematics and lensing, the average of the 
estimator is given as 



(oebCL)) 



a(L) + 

lu(L)+A eb (L) 



(27T) 



2 SebFeb<KL), 
(Al) 



where the first term on the right hand side is the desired 
rotation field, the second term is the estimator bias from 
instrumental systematics, and the third term represents 
the bias from lensing. The variance of the estimator can 
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be written as 

= A EB {L)A EB {L') x 

= (2tt) 2 (5 d (L + L') x 

C aa (L) + C""{L) + N EB , EB (L) + N% BtEB {L) 



+N%% EB (L)+N^% B (L) + 



(A2) 



where L = li + I2, C aa (L) is the cosmological rotation 
power spectrum, and C" U {L) is the rotation systematics 
power spectrum. The terms N EBtEB (L), N EB EB (L), 

N EB EB (L), and N EBEB (L) arc the estimator Gaus- 
sian noise, the first order non-Gaussian estimator noise, 
the first order systematics noise of instrumental rotation, 
and the first order lensing induced non-Gaussian noise, 
respectively. Like the lensing quadratic estimator, the 
Gaussian noise comes from the disconnected part of the 
four-point function, while non-Gaussian noise N^ EB EB , 



N 



EB,EB 



, and N 



(le 

EB.EB 



comes from the connected part. 



We note that the Gaussian noise term also includes ro- 
tation systematic effects implicitly since instrumental 
systematics bias the measured rotation power spectrum 
with C UUJ (L) with as we have shown, note that we as- 
sume no cross correlation term C au} {L). The ellipses in 
Eq. (|A2[) stands for higher order terms. The first order 
non-Gaussian noise can be written as 



iV 



EB.EB 



(L) = 



A 2 EB (L) 



d 2 h 



(2tt) 2 J (2tt) 2 ' 

c^ E ol[ E lo^ hl w§(h,-h)w§(i' 2 ,-i[) 



+C^y 2 \W B c (h,-l> 1 )W§(l> 2 . 



-li) 



(A3) 



where X = {a,u>,(j> len }, the window functions 
W^(li,l a ) = W^(li,la) = 2cos[2(^ 2 - <p h )\, and 
W|' e "(li,l 2 ) = sin2(( / 3 ll - ip h ){L ■ h). We use this 
equation to numerically compute the systematic-induced 
estimator noise for the rotation systematics. Among 
these extra covariance noise, A^°g EB , and N E L B EB 
are always smaller than the estimator Gaussian noise 
N EB , EB {L), and N EB EB can be in some cases compa- 
rable to N EB:EB (L). We use Fig. [5] to illustrate when 
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Neb eb S oes to tnc same level as Neb,eb(L) under cer- tain experiment configuration. 



